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Alerting Abstract: WO AI 

NOVELTY - The position of the imaging source (7) is found fiom a reference repository, and a predefined model matching the mean form of the object selected . This is 
repeated until variation between the projection of contours of the two-dimensional image and the model surface are minimal. Model orientation and position are selected from 

a repository, and a three-dimensional deformation selected. 

USE - Reconstruction of x-ray images from partial data to visualize bones, for use in medicine or archaeology. 

ADVANTAGE - Allows construction of the image of the bone surfaces, without cost of using an X-ray scaimer and requiring only low levels of radiation exposure for 
patients. 

DESCRIPTION OF DRAWINGS - The drawing shows a diagram of projections fiom the object 
7Imaging source 

1 Object 

Main Drawing Sheet(s) or Clipped Structure(s) 
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RECONSTITUTION STATISTIQUE DE SURFACES EN TROIS DIMENSIONS** 

Assignee: UNIVERSITE JOSEPH FOURIER, 621, avenue Centrale, Boite postale 53. 38041 Grenoble Cedex 9, FR 

Inventor: FLEUTE, Markus, 6, rue des Peupliers, F-38400 Saint Martin d'Heres. FR LA\'ALLEE. Stephana, 4, rue Augereau, F-38000 Grenoble, FR DESBAT, Laurent, 1, 
rue Ainedee Morel, F-38000 Grenoble, FR 

Agent: de Beaumont, Michel, Cabinet Michel de Beaumont, 1 , rue Champollion. 3 SOOO Grenoble, FR 
Language: FR 

Application: EP 2000962642 A 20000914 (Local application) WO 2000FR2546 A 20000914 (PCX Application) 
Priority: FR 199911848 A 19990917 
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Original Abstract: The invention concerns a method and a system tor three-dimensional reconstmction of an image representing the surface contours of at least an object (1), 
from at least a two-dimensional view of said object obtained by X-ray, which consists in: determining the position of the photographing source (7) in a reference repository; 
selecting a predefined model constituting a mean fonn of the object, and repeating the process until the contours of the model are such that the variations between the 
overhead projection rays of the contours of the two-dimensional image fiom the source and the model surface are minimal; selecting an orientation and a position for the 
model in the reference repository; then selecting a deformation of the model to modify its contours in three dimensions. jEP 1222636 Bl (Update 200421 E) 
Publication Date: 200403 1 7 

**DREIDIMENSIONALE STATISTISCHE FLACHENREKONSTRUKTIGX THREE-DIMENSIONAL STATISTIC RECONSTRUCTION OF SURFACES 
RECONSTITUTION STATISTIQUE DE SURFACES EN TROIS DIMENSIONS** 

Assignee: UNIVERSITE JOSEPH FOURIER, 621, avenue Centrale, Boite postale 53, 38041 Grenoble Cedex 9, FR 

Inventor: FLEUTE, Markus, 6, rue des Peupliers, F-38400 Saint Martin d'Heies, FR LAVALLEE, Stephane, 4, rue Augereau, F-38000 Grenoble, FR DESBAT, Laurent, 1, 
rue Amedee Morel, F-38000 Grenoble, FR 

Agent: de Beaumont, Michel, Cabinet Michel de Beaumont, 1 , rue Champollion, 3 8000 Grenoble, FR 
Language: FR 

Application: EP 2000962642 A 20000914 (Local application) WO 2000FR2546 A 20000914 (PCT Application) 

Priority: FR 19991 1848 A 19990917 

Related PubUcation: WO 200 1 022368 A (Based on OPI patent ) 
Designated States: (Regional Original) DE FR 
Original IPC: G06T-1 1/00(A) 

Current IPC: A61B-19/00(R,I,M,EP,20060101,20051008A) A61B-19/00(R,IAI,EP,20060101,20051008,C) 
Current ECLA class: A61B-19/00N A61B-19/00N6T 

Current ECLA ICO class: K61B-19:00E2P K61B-19:00N10 K61B-19:00N10C K61B-19:00N6T K61B-19:00N6T1L K61B-19:00N6T4 K61B-19:00N6T5 K61B-19:00U6 
Claim: 1 .Verfahren zur Rekonstitution eines die Oberflachenkonturen wenigstens e ines Objekts (1) wiedergebenden Bildes in drei Dimensionen, ausgehend v on wenigstens 

einer mit Rontgenstrahlung gemachten zweidunensionalen Au fnahme dieses Objekts. **dadurch gekennzeichnet, dass** das Verfahren umfasst: * Bestunmen der Lage bzw. 
Stelltmg der Aufnahmestrahlungsquelle (7) in einem Bezugs-Referentialsystem: * Wahlcn \\ cnig.-.icns eines statistischen Modells, welches eine durchsc hnittliche mittlere 
Form des Objekts und seine Hauptdeformationen r elativ bezuglich dieser mittleren diiich.schnittlichcn Form definiert , wobei das statistische Modell ausgehend von einer 
Population von Objekten dcsselben Typs bestimmt wird. fur welche man die samllichc n Objcktcn gcmcinsame statistische Kon'e.spondenz-Entsprechung sucht ; sowie * 
iterative Duiohfiihnmg der folgenden Schritte, bis die Konturen sol cher Art sind, dass die Abstande zwischen den Ruckprojektionsstrahl en der Konturen des 
zweidimensionalen Bildes von der Strahlungsquel le und der Modelloberflache minimal sind, zur Erlangung einer Korre spondenzentspiechung zwischen dem Modell und 
dem Bild: * Wahlen einer Orientierung bzw. Ausrichtung und einer Stellung bzw . Lage des ModeUs in dem Bezugsreferentialsystem, indem man das statistische Modell 
einer stanen Transformation unter Anderung s einer Stellung/Lage und/oder seiner Orientierung /Ausrichtung unt erwirfl, sodann * Wahlen einer Deformation des Modells, 
indem man das statistische Modell einer nicht-starren Deformation unter Veranderung seiner O berflachenkonturen unterwirfl. 1 .A method for restoring a three-dimensional 
image representing the surfa oe contours of at least one object (1), based on at least one two-dimen sional X-ray view of this object, **characterized in that** it consist s of: * 
determining the position of the shooting source (7) in a reference re ferential system; * selecting at least one statistical model defming an average shape (2 1) of the object and 
its main deformations with respect to this avera ge shape, the statistical model being calculated fiom an object popul ation of a same type for which the statistical 
correspondence common to all objects is searched; anditeratively, until the contours of th e model are such that the inter\'als between back -projection rays ol~l he image 
contours in two dunensions from the source and the model sur face are minimum, to obtain a correspondence between the model and th e unage: * selecting an orientation and 
a position of the model in the reference referential system by submitting the statistical model, successively , to a rigid transformation modifying its position and/or its orienta 
tion, then * selecting a deformation of the model to modify its contours in three dimensions, by submitting the statistical model to a non-rigid deform ation modifying its 
surface contours. l.Procede de reconstitution d'une image en ttois dimensions representant les contours surfaciques d'aumoins unobjet (1), apartird'aumoinsune vue en 
deux dimensions de cet objet prise aux rayons X, **caracteri se en ce qu'**il consiste a: * determiner la position de la source (7) de prise de vues dans un refe rentiel de 
reference; * selectionner au morns un modele statistique defrnissant une forme moy enne de I'objet et ses deformations principales par rapport a cette f orme moyenne. le 
niodele slalislique etant calcule apartir d'une popu lation d'objets de meme type pour laquelle on recherche la cortespond ance stalisliquc coninuine a lous k\s ohjets: el de 
laeon ileralive juscju'a ce cjue les contours du modele soient tels que les ecarts, ent re des rayons de retroprojection des contoru's de I'unage en deux dime nsions depuis la 
source et la surface du modele, soient minunaux, pou r obtenir une correspondance entre le modele et I'unage: * selectionner une orientation et une position du modele dans le 
refere ntiel de reference en faisant subir au modele statistique une transfo rmation rigide modifiant sa position et/ou son orientation, puis * selectionner une deformation du 
modele en faisant subir au modele sta tistique une deformation non rigide modifiant ses contours surfacique s. 
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Related Publication: WO 2001022368 A (Based on OPI patent ) 
Original IPC: G06K-9/00(A) 
Cmrent IPC: G06K-9/00(A) 

Cmrent ECLA class: A61B-19/00N A61B-19/00N6T 

Cmrent ECLA ICO class: K61B-19:00E2P K61B-19:00N10 K61B-19:00N10C K61B-19:00N6TK61B-19:00N6T1LK61B-19:00N6T4 K61B-19:00N6T5 K61B-19:00U6 

Cmrent US Class (main): 382-132000 

Cmrent US Class (secondary): 382-128000 382-154000 

Original US Class (main): 382132 

Original US Class (secondary): 382128 382154 

Original Abstract: The invention concerns a method and a system for three-dimensional reconstruction of an image representing the smface contoms of at least an object 
(** 1 **), from at least a two-dunensional view of said object obtained by X-ray, which consists in: determining the position of the photographing somce (**7**) in a reference 
repository; selecting a predefined model constituting a mean form of the object, and repeating the process until the contours of the model are such that the variations between 
the overhead projection rays of the contoms of the two-dimensional linage from the source and the model surface are minimal; selecting an orientation and a position for the 
model in the reference repository; then selecting a deformation of the model to modify its contours in three dunensions. 

Claun: The invention claimed is: 1 . 1 . A method for restoring a three-dimensional image representing the su rface contours of at least one object, based on at least one two- 
dunens ional X-ray projection of this object, the method comprising the steps of: * determining the position of the shooting source in a reference refere ntial system; * 
selecting at least one statistical model defLning an average sh^e of the object and its main deformations with respect to this average sh ^e, the statistical model being 
calculated from an object population of a same type for which the statistical correspondence common to al 1 objects is searched; and * selecting an orientation and a position 
of the model in the reference referential system by submitting the statistical model, successively , to a rigid transformation modifying its position and/or its orienta tion, then * 
selecting a deformation of the model to modify its contom's in three dimensions, by submitting the statistical model to a non-rigid deform ation modifying its smface contours, 
* wherein the orientation and deformation selecting steps are performed iteratively, imtil the contours of the model are such that the inter vals between back-projection rays of 
the image contours in two dimens ions from the source and the model surface are minimimi, in order to o btain a correspondence between the model and the image. 

WIPO 

Publication Nmnber: WO 2001022368 Al (Update 200158 B) 
Publication Date: 20010329 
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Augereau, P-38000 Grenoble. PR Residence: PR Nationality: PR -(only US )- nF.SH,.'iT, Laurent, 1 , rue Amedee Morel, P-38000 Grenoble, PR Residence: PR Nationality: 
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Grenoble. I R Residence: I R Nationality: PR DESBAT, Laurent, 1, rue Amedee Morel, F-38000 Grenoble, FR Residence: FR Nationality: FR 
Agent: DE BEAUMON T, Michel, 1, rue ChampoUion, F-38000 Grenoble, FR 
Language: FR (29 pages, 5 drawings) 

Application: WO 2000FR2546 A 20000914 (Local application) 
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Original Abstract: The invention concerns a method and a system for three-dimensional reconstruction of an image representing the surface contours of at least an object (1), 
from at least a two-dimensional view of said object obtained by X-ray, which consists in: determining the position of the photographing soince (7) in a reference repository; 
selecting a predefined model constituting a mean form of the object, and repeating the process imtil the contours of the model are such that the variations between the 
overhead projection rays of the contoms of the two-dimensional image from the source and the model surface are minimal; selecting an orientation and a position for the 
mode I i n t he reference repository; then selecting a deformation of the model to modify its contom's in three dimensions. L'invention conceme \m precede et \m systeme de 
reconstitution d'une unage en trois dimensions representant les contoms smfaciques d'aumoins im objet (1). a partie d'au moins une vue en deux dimensions de cet objet prise 
aux rayons X, dans lesquels on determine la position de la som'ce (7) de prise de vues dans im referentiel de reference, on selectionne un modele predetlni constituant ime 
forme moyenne de I'objet et, de facon iterative jusqu'a ce que les contours du modele soient tels que les ecarts entre des rayons de retroprojection des contoms de I'image en 
deux dimensions deptus k source et ia surface du modele soient minimaux, on selectionne tme orientation et tme position du modele dans le referentiel de reference, ptiis on 
selectionne une deformation du modele pour modifier ses contours en trois dunensions. 
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RECONSTITUTION DE SURFACES EN TROIS DIIUENSIONS PAR UTILISATION DE MODELES STATISTIQUES. 

_ L'inventlon concerne un procede et un systeme de re- 
constitutlon d'une image en trois dimensions representant 
les contours surfaciques d'au moins un objet (1), a partir 
d'au moins una vua an daux dimensions da cat objet prise 

aux rayons X, dans lesquels on determiner la position de la 
source (7) de prise de vues dans un referentiel de referen- 
ce, on selectionne un modele predefini constituant une for- 
me moyenne de I'objet et, de fa^on iterative jusqu'a ce que 
les contours du modele soient tels que les ecarts entre des 
rayons de retroprojection des contours de I'image en deux 
dimensions depuis la source et la surface du modele soient 
minimaux, on selectionne une orientation et une position du 
modele dans le referentiel de reference, puis on selection- 
ner une deformation du modele pour modifier ses contours 
en trois dimensions. 
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RECGNSTITUTION DE SURFACEiS EST TROIS DIMENSIONS FAR XITILISATIQN DE 

MODELES STATISTIQDES 

La presents invention conceme la reconstitution d'itta- 
ges en trois diinensions et, plus particuli^retnent, la reconstitu- 
tion des contours surfaciques d'une image k partir de vues, mitne 
incompletes, en deux ditnensions . 
5 Un exenple d ' application de la pr6sente invention est 

la reconstitution d' images repr^sentant des elements osseux i 
partir d' images en deux ditnensions prises, par exeirple, aux 
rayons X. De telles images peuvent servir, par exettple, a simuler 
una intervention chirxirgicale en permettant au praticien de pr6- 

10 visualiser les positions respectives des os au niveau d'une arti- 
culation, par exetiple, en cas de plastie ligaitientaire ou de mise 
en place de prothese. 

Aujourd'hui, pour permettre au praticien de visualiser 
une image en trois dimensions a partir d' images en deux dimen- 

15 sions, on est contraint d'utiliser un scanner t rayons X effec- 
tuant des tomographies. A partir de ces tomographies, on sait 
reconstituer une image en trois dimensions. Soit il s'agit d'une 
technologie dite de tomographies en trois dimensions dans 
laquelle un grand norribre de vues (de I'ordre de 200) sont prises 

20 au moyen d'vm scanner a rayons X selon diff^rentes positions, et 
les contours de I'os sont alors reconstituis par traiteraent 
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d' image partir de ces vues au scanner. Soit il s'agit d'me 
technologie dite de tomographies en deux dimensions dans laquelle 
vox grand norribre de coupes sont prises au moyen d'lxn scanner 
rayons X transversalement a I'os. On peut alors reconstituer 
5 1' allure et la structure de I'os. 

La technique du scanner donne de bons r^sultats mais 
est d'une mise en oeuvre lourde et couteuse. En effet, 1' utilisa- 
tion d'm scanner permet d'obtenir un ensemble d' images en deux 
dimensions foumissant non seulement des informations sur le 

10 contour mais ggalement sur I'interieur de I'os. Or, dans de ncm- 
breuses applications, seule la connaissance du contour surfacique 
de I'os ou de I'objet est necessaire. 

Un autre exenple d' application de la presente invention 
est la reconstitution d'ossements incottplets, par exenple, en 

15 arch6ologie. \Jne reconstitution d' images en trois dimensions peut 
permettre de retrouver de fagon quasi -parfaite la forme origi- 
nelle de I'os rairae si celui-ci est decouvert de fagon incoirpldte. 
Dans une telle application, des prdbl&nes similaires k ceux expo- 
sis ci-dessus en relation avec la simulation d ' interventions chi- 

2 0 rurgicales se posent. En particulier, il est souvent utile de 

connaitre la forme de I'os sans avoir i se prioccuper de sa 
structure interne. 

Un. autre inconvenient des techniques connues est 
qu'elles imposent une dose de radiations importante pour le 
25 patient ce qui n'est pas souhaitable. Si cet inconv&iient est 
moins sensible dans le doraaine de 1 ' archiologie ou c'est le cout 
qui pr&Jomine, il est particuli^rement ggnant dans la simulation 
d' interventions chirurgicales . 

La presente invention s 'applique plus particulierement 

3 0 ^ la reconstitution d' images se rapportant k des objets identi- 

fies, c'est-a-dire, dont on connait a I'avance la forme generale. 
Par exeirple, pour vin os, il faut au prealable decider de quel os 
il s'agit. 

La presente invention vise a proposer un nouveau pro- 
35 cede de reconstitution d' images en trois dimensions qui pallie 
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les inccaiv6nients des techniques connues. L' invention vise, en 
particulier, S. proposer una solution qui ne nicessite pas I'em- 
ploi couteux d'un scanner i. rayons X. 

La pr§sente invention vise egalement a proposer me 
5 solution qui soit connpatible avec me exposition tninimale aux 
rayons X ou Equivalents. 

L' invention vise en outre a minimiser le nombre de vues 
en deux dimensions necessaires pour reconstituer 1' image en trois 
dimensions . 

10 Pour atteindre ces dbjets, la presente invention pre- 

voit m precede de reconstitution d'me image en trois dimensions 
representant les contours surfaciques d'au moins m objet, ^ par- 
tir d'au moins me vue en deux dimensions de cat objet prise aux 
rayons X, caract^rise en ce qu'il consist e 4 : 

15 determiner la position de la source de prise de vues 

dans m rifirentiel de reference ; 

sElectiomer m module prgdefini constituant me forme 
moyeme de I'dbjet / et 

de fagon iterative jusqu'S. ce que les contours du 

20 module soient tels que les hearts, entre des rayons de r6tropro- 
jection des contours de 1' image en deux dimensions depuis la 
source et la svirface du modile, soient minimaux : 

sElectionner me orientation et me position du module 
dans le rifirentiel de reference, puis 

25 s^lectiomier me deformation du raoddle pour modifier 

ses contours en trois dimensions. 

Selon m mode de realisation de la presente invention, 
le moddle est obtenu a partir d'me population d ' objet s pour 
laquelle on recherche la correspondance statistique comrame ^ 

3 0 tous les objets pour determiner me forme moyeme et les deforma- 
tions principales par rapport a cette forme moyeme, de fagon a 
disposer d'au moins m modele statistique. 

Selon m mode de realisation de la presente invention, 
les etapes de selection iteratives consistent a faire siibir au 

35 modele statistique, successivement, me transformation rigide 
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itradifiant sa position et/ou son orientation et une deformation 
non-rigide tnodifiant ses contours surfaciques. 

Selon un mode de realisation de la presents invention, 
les contours de 1 ' image en deux dimensions sont obtenus automat i- 
5 queraent en projetant le modele dans le plan de 1' image en deux 
dimensions, et en diformant les contours projetes de fagon a les 
faire coincider avec les points de fort gradient en niveaux de 
gris de 1' image en deux dimensions. 

Selon vn mode de realisation de la presents invention, 
10 la determination automat ique des contours de 1' image en deux 
dimensions est effectuee de fagon iterative, chaque iteration 
etant intercalee entre deux iterations successives des etapes de 
selection. 

Selon un mode de realisation de la presente invention, 
15 on determine, dans le referent iel de reference, des coordonnees 
en trois dimensions de points de I'dbjet, de fagon k disposer de 
points de reference suppiementaires pour les etapes de selection 
iteratives de position, orientation, et deformation. 

Selon im mode de realisation de la presente invention, 
20 on utilise plusieurs images en deux dimensions pour lesquelles 
les positions respectives de la source de prise de vues sont tou- 
tes determinees dans le referent iel de reference, et on effectue 
les etapes de selections iteratives en tenant cottpte des rayons 
de retrqprojection des contours de toutes les images en deux 
2 5 dimensions . 

Selon un mode de realisation de la presente invention, 
le nortibre d' images utilise est fonction de la precision souhai- 

tee. 

Selon un mode de realisation de la presente invention, 
30 la surface du modele est constituee d' elements de triangle, les- 
dits ecarts etant mesures par rapport a des points de certaines 
aretes constituant des generateurs du contour en trois dimen- 
sions . 

Selon un mode de realisation de la presente invention, 
35 le precede est applique a la reconstitution des contours surfaci- 
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ques de plusieurs objets li€s entre evtx par des relations de 
transfortnaticn rigide et/ou elastique. 

Selon un mode de realisation de la presente invention, 
le precede est applique a la reconstitution d' images d'os. 
5 L' invention conceme ^galement \in systeme de traitement 

d' images, cortportant des moyens pour la mise en oeuvre du proced§ 
de reconstitution d' images en trois dimensions. 

Ces objets, caracteristiques et avantages, ainsi que 
d'autres de la presente invention seront exposes en detail dans 
10 la description suivante de modes de realisation particuliers 
faite a titre non-limitatif en relation avec les figures jointes 
parmi lesquelles : 

la figure 1 repr§sente un mode de realisation d'un sys- 
tdme de reconstitution d' images en trois dimensions i. partir d'un 
15 dispositif de prise de vues en 6evx dimensions selon la pr§sente 
invention ; 

la figure 2A illustre la prise d'une vue en deux dimen- 
sions pour la mise en oeuvre du proc§d6 de reconstitution d' ima- 
ges en trois dimensions selon la presente invention ,- 

2 0 la figure 2B illustre la prise de trois vues en deux 

dimensions pour la mise en oeuvre du proc6d6 de reconstitution 
selon la present e invention / 

les figures 3A, 3B et 3C illustrent trois etapes carac- 
teristiques du precede selon m nuDde de realisation de la prg- 
25 sente invention / 

la figure 4 illustre un mode de determination de g&ie- 
rateurs de contour selon la presente invention ; et 

la figure 5 illustre vn mode de determination pref^re 
selon 1' invention de I'^cart entre un module en trois dimensions 

3 0 et une vue en deux dimensions. 

Les memes elements ont ete designes par les mSmes refe- 
rences aux differentes figures qui ne sont pas a I'echelle. Pour 
des raisons de clarte, seuls les elements qui sont necessaires a 
la cotrprehension de 1' invention ont ete reprgsentes aux figures 
35 et seront decrits par la suite. En particulier, les moyens de 
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traitement d' images cortprenant les calculateurs informatiques 
n'ont pas ete repr§sentes et ne seront pas decrits de fagion 
detainee, car leur realisation est a la portee de I'honme du 
Tn§tier a partir des indications fonctionnelles donnees ci-apres. 
5 Une caracteristique de la presente invention est de 

rechercher la position et 1 'orientation d'un objet a reconsti- 
tuer, et dont on connait au moins une vue en deux dimensions, a 
partir d'une base de donnees contenant des modeles de cet objet. 
Quand plusieurs vues en deux dimensions sont utilisees, ces vues 

10 sont toutes r^ferencees dans un meme referentiel. Ainsi, 1' inven- 
tion conceme la reconstitution d'une image reprSsentant \m objet 
d#j§. identifi^ et dont on peut disposer de modeles de taille 
et/ou de forme differentes. 

Une caracteristique d'un mode de realisation pr^f^r^ de 

15 la presente invention est d'utiliser au moins un moddle statisti- 
que d§forraable, 6tabli k partir de la base de donnees, pour 
reconstituer la forme en trois dimensions de 1' objet. Ainsi, 
1' invention pr^voit de dSfinir, avant toute reconstitution, une 
base de donnees contenant des modules en trois dimensions de 

20 1' objet k reconstituer ou, de preference, un ou plusieurs modeles 
statistiques dgformables i partir de cette base de donnees. 

La figure 1 reprisente une vue schiraatique d'xin systeme 
de reconstitution d' images en trois dimensions selon un mode de 
realisation de la prisente invention. L'exetiple de la figure 1 

25 conceme la reconstitution de 1' image d'un os 1 k partir de 
radiographies en deux dimensions. En figure 1, on a represents 
schematiquement et partiellement le corps d'lan patient p dont la 
jattibe cent lent I'os l que I'on souhaite visualiser. La jambe du 
patient p (done I'os 1) est, par exetrple, pose sur une table 2. 

3 0 Une potence 3 porte un dispositif 4 de localisation en trois 
dimensions qui peut etre optique, magnetique, mecanique ou ultra- 
sonore et qui repere la position et 1 ' orientation de multiples 
capteurs-emetteurs constitues, par exenple, de diodes infra- 
rouges, de reflecteurs, d'emetteurs magnetiques, a ultrasons, 

35 etc. Un tel localisateur est parfaitement connu et ne sera pas 
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detains plus avant. On notera siTrplement que, pour garantir une 
localisation conrecte du patient (de I'os 1), celui-ci est gene- 
ralement 6galenient Equips d'un capteur-emetteur 18 detectable par 
le localisateur 4. En effet, selon 1' invention, toutes les vues 
5 doivent pouvoir etre eacploitees dans un mStne referentiel (de 
r§f §rence) , qui est associe a I'objet que I'on reconstitue. 

Plusieurs capteurs-emetteurs peuvent etre fixes sur le 
systeme de radiologie, proche de la source 7 (capteur-emetteur 
19) ou proche d'un detecteur d' image 9 (capteur-emetteur 5), de 

10 fagon a reperer la position du systeme de radiologie par rapport 
au referentiel de reference du capteur-emetteur 18. 

Dans certains cas, les capteurs-emetteurs 5 et 19 sont 
difficileraent reperables en raison de leur eloignement ou de la 
presence d'objets parasites dans le champ de ttiesure du localisa- 

15 teur 4, Dans ce cas, un capteur-emetteur 20 peut it re ins t alii 
sur la table 2 dans le chaitp de mesure du localisateur 4. On 
positionne alors le systdme radiologique dans de bonnes condi- 
tions de raesure et on rep^re les positions des capteur-^metteurs 
5 et 19 par rapport au capteur-emetteur 20 vine fois pour toutes 

20 (on ne renouvelle cette ^tape que si on bouge le systdme de 
radiologie dans son ensettible) . 

Le syst&tie radiologique etant ^quipi de codeurs angu- 
laires sur ses axes, comme cela sera dStailie plus loin, les 
changements de position relative du syst&tie radiologique sont 

25 raesur^s ^ I'aide de ces codeurs et peuvent ainsi etre report6s 
dans le referentiel du capteur-emetteur 20. Pour chaque prise de 
vue, on Ttiesure la relation geometrique entre les capteurs- 
emetteurs 20 et 18 et, par ce biais, on reporte la position du 
systeme radiologique dans le referentiel de reference du capteur- 

3 0 emetteur 18. Ainsi, toutes les radiographies sont calibrees dans 
le referentiel du capteur-emetteur 18 qui peut etre mobile. 

Le dispositif 6 est, par exetrple, constitue d'une 
source 7 a rayons X portee par une premiere extremite d'un bras 8 
en demi-cercle dont 1' autre extremite est destinee a recevoir le 

3 5 film 9 d ' impression de la radiographie, ou un capteur electron!- 
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que Equivalent tel qu'un aitplificateur de brillance, ou un d6tec- 
teur plat au silicium amorphe. La relation entre le bras 8 et la 
table 2 est telle que cette demidre se trouve entre la source 7 
et le capteur 9. Le bras 8 est tnonte a rotation autour d'un axe 
5 10, motorise ou a deplacement tnanuel, et porte par me potence 
11. L' ensemble de prise de vues en deux dimensions peut ainsi 
toumer autour de I'os 1 pour effectuer le nombre desire de 
radiographics de celui-ci. Le cas echeant, le localisateur opti- 
que 4 peut etre associe a un dispositif de codage angulaire de la 

10 position de I'axe 10. 

L' ensemble est pilote par vn systeme informatique, par 
exerrple, un ordinateur 12 associe, de preference, a un ecran de 
visualisation 13 . A la figure 1 , on a symbolise par des liaisons 
unifilaires 14, 15, 16 et 17, les bus infontiatiques d'Echange de 

15 signaux Electriques de commande et de donnees entre 1' ordinateur 
12 et, respectivement, le localisateur 4, la source 7, le inoteur 
ou le codeiiT cptionnel de I'axe 10 et le capteur 9. 

Les figures 2A et 2B illustrent deux modes de prises de 
vues selon deux modes de realisation de la pr^sente inventioa. 

20 Aux figures 2A et 2B, les vues en deux dimensions ont §t6 repri- 
sent6es par les plans respectifs dans lesquelles elles sont pri- 
ses, c'est-^-dire par 1' allure de la svurface du capteur 9 (figure 
1) lors des prises. Cela correspond aux images radiographiques en 
deux dimensions ricup^r^es par le syst§me informatique. 

25 La figure 2A repr6sente le cas d'une seule prise de vue 

dans un plan PI, foumissant une image II de I'os 1. La position 
de la source 7 a ete syttibolisEe par un point ou convergent des 
rayons r de r^trqprojection des quatre coins du plan PI. 

La figure 2B repr^sente le cas d'une triple prise de 

30 vues dans des plans PI, P2 et P3, foumissant trois images II, 12 
et 13 en deux dimensions de I'os 1. En figure 2B, la position de 
la source n'est pas constante, celle-ci est differente pour cha- 
que prise de vue. Toutefois, grdce au localisateur 4, toutes les 
positions de la source sont connues dans le referentiel de refe- 
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rence. Par soucis de clart6, seuls les rayons rl des quatre coins 
du plan PI dans la position de I'itnage II ont ete representes. 

On notera que la source de prises de vues peut subir, 
entre devix itrages en deux dimensions, d'autres mouvements que 
5 dans m rngme plan comme illustre par les figures 1 et 2B. En 
d'autres temes, le bras 8 du systeme de prises de vues 6 peut 
presenter plus de deux degres de liberte, dent chacun pourra etre 
muni d'un dispositif de codage angulaire. Selon un mode de reali- 
sation prefere, le bras 8 est monte sur deux axes de rotations 

10 horizontal et vertical et un axe de translation vertical. 

On notera egalement que les images en deux dimensions 
peuvent n'etre que partielles. Par exeitple, dans 1 ' application 
aux radiographies , celles-ci peuvent etre interpretees par I'ope- 
rateur pour valider les contours des vues en deux dimensions IL 

15 prendre en corrpte pour la reconstitution. Une telle interpreta- 
tioa n'est pas g§nante en raison du faible nottibre de vues en deux 
dimensions n§cessaire selon 1' invention (moins de dix en g§n^- 
ral) . La validation des contours dans le systeme infomiatique 
peut §tre effectu^e, par exeraple, au raoyen d'lane souris, d'm 

2 0 crayon cptique, d'm 6cran tactile, ou Equivalent, de fa?on clas- 
sique pour un enregistrement de contour sur une image en niveaux 
de gris. 

Selon un autre mode de realisation prSfgrg, la determi- 
nation des contours sur les images radiographiques est autcmati- 

25 s^e en raettant en oeuvre un procSde dit de recalage en deux 
dimensions. Un tel proc6d6 consiste k determiner le contour de 
fagon automatique par analyse des gradients les plus forts en 
niveaux de gris. Ce procede est decrit, par exemple, dans la 
these de Gelu lonescu, presentee publiquement le 4 deceiribre 1998 

30 a 1' University Joseph Fourier de Grenoble (France), et ayant pour 
titre "Segmentation et recalage d' images echographiques par uti- 
lisation de connaissances physiologiques et morphologiques" . 
Selon 1' invention, on met en oeuvre ce procede de recalage en 
deux dimensions en combinaison avec une projection, sur 1' image 

35 en deux dimensions, du modele statist ique deforraable de I'objet 
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dont on souhaite reccnstituer les contours surfaciques en trois 
dimensions. On effectue una analyse iterative d'lon contour pro- 
jet e du tnodele sur 1' image en deux dimensions jusqu'a obtenir, 
apres plusieurs passes (par exenple de 3 a 10) , une identite 
5 entre le contour du modele projetS qui sera decrit ci -apres et 
dont 1 ' orientation et la deformation auront ete adaptes par le 
procede de 1' invention, et le contour determine par analyse des 
gradients les plus forts en niveaux de gris. Cette cotribinaison 
sera mieux cottprise par la suite en relation avec la figure 5. 

10 Pour sittplifier la presente description, on considerera 

ci-apres que les contours des vues en deux dimensions sont 
connus. En pratique, les etapes du procede de 1' invention seront 
effectuees en boucle, y coopris le recalage en deiox dimensions, 
jusqu'a obtenir un resultat satisf aisant . 

15 Selon 1' invention, une fois que la ou les vues en deux 

dimensions ant 6t§ dbtenues, il reste §l diterminer, par recherche 
de correspondance dans la base de donnges, la forme et la taille 
de I'objet (ici, I'os) en trois dimensions. 

Pour ce faire, 1' invention prgvoit de rechercher le 

2 0 modile pour lesquelles les distances separant les contours de 

chaque vue en deux dimensions et la surface du module sont mini- 
males. 

Dans le mode de rialisation pr#f6ri ou on utilise un 
modele statistique d^formable (definissant ainsi une famille de 
25 modules) de I'dbjet i reconstituer, on recherche ^ s'approcher le 
plus possible, par deformation iterative de ce module, de la 
forme dont les contours sont tels que, corame on le verra par la 
suite en relation avec la figure 5, les rayons de retroproj ection 
s'appuyant sur les points des contours de 1' image (ou des images) 

3 0 en deux dimensions depuis la source (ou depuis leurs sources 

respectives) sont tous tangentiels a la surface du modele, etant 
entendu qae les differents contours en deux dimensions sont 
reperables dans le referentiel de reference. 

La realisation du modele statistique deforraable en lui 
35 meme ne fait pas I'objet de la presente invention et est parfai- 
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tement classique. On se bomera t rappeler que la realisation 
d'un tel nradile fait g&ieralement appel a une recherche d'une 
forme itioyenne d'une population d'objets de meme type (par exam- 
ple, des fimurs) const ituant la base de donnees, suivie d'une 
5 analyse en composantes principales pour determiner les d^fontia- 
tions principales (essentielles) a appliquer a cette forme 
mpyenne et dbtenir ainsi le modele statistique. 

Une tti^thode de determination d'un modele statistique a 
partir d'une famille d'echantillons est decrite, par exettple, 

10 dans 1" article "Building a Con^plete Surface Model from Sparse 
Data Using Statistical Shape Models: Application to Conputer 
Assisted Knee Surgery" de Markus Fleute et Stephane Lavallee, 
paru dans medical image computing and computer-assisted 
INTERVENTION - MICCAI'98, Springer -Ver lag LNCS Series, pages 880- 

15 887, en Octobre 1998, et dont le contenu est incorpori ici par 
reference. On notera que la precision de la reconstruction des 
contours surfaciques qp6r§e par 1' invention depend des 4chantil- 
lons utilises pour construire la base de donnees du module sta- 
tistique. Par exeraple, si vm module statistique de f6raur est crit€ 

20 partir d'lme population d'ichantillons ayant des formes norma- 
les (sans pathologie) , on pourra reconstituer avec precision des 
formes normales mais la precision sera limit^e si le f6mur radio- 
graphi§ t me forme pathologique . Par centre, si le module stati- 
stique est cr€§ k partir d'une grande population d'echantillons 

25 contenant k la fois des formes normales et des formes pathologi- 
ques, on pourra reconstituer, avec precision, des objets normaux 
et pathologiques . 

Selon 1' invention, une fois que I'on dispose des images 
en deux dimensions, on recherche le cas icheant k quel modele 

30 contenu dans la base de donnees ces images se rapprochent le plus 
(par exenple, si plusieurs types d'os ou si plusieurs modeles 
statistiques d'm meme os sont disponibles) , puis on applique des 
transformations a ce modele pour s'approcher le plus possible 
d'une image en trois dimensions, telle que les rayons de retro- 

35 projection des images jusqu'a leurs sources respect ives sont tan- 
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gentiels aux surfaces de 1' image en trois dimensions, plac6e sur 
le trajet de ces rayons et dont ils sont respectivement les plus 
proches . 

De pr§f France, on commence par faire subir au modele 
5 initial au repos, c'est-a-dire correspondant ^ la forme moyenne, 
une transformation dite rigide, c'est-a-dire, ne touchant qu'a 
son orientation et sa translation dans I'espace. Puis, quand on a 
obtenu la position pour laquelle les rayons de retroprojection 
sont tous §. une distance minimale des contours du modele initial, 

10 on fait subir une deformation non-rigide au modele, c'est-^-dire 
que sans modifier son orientation, on modifie sa forme a partir 
des donnees contenues dans la base statistique en modifiant les 
coefficients des modes principaux du modele, jusqu'a obtenir la 
forme dont les contours se rapprochent le plus des rayons de 

15 retroprojection. Le cas 6ch6ant, le modile est prepositionn^ de 
fagon tr§s approximative par I'cp^rateur par rapport aux vues en 
deux dimensions afficMes $l I'ecran. On notera que les moyens 
(informatiques) de traitement d' images utilises sont classiques 
dans leur structure et n'ont done pas besoin d'itre ditaillis. 

20 Les figures 3A h 3C illustrent ces deux Stapes du pro- 

c^6 de 1' invention. La figure 3A reprisente une vue en trois 
dimensions d'un modele 21 avant toute deformation. II s'agit, par 
exemple, d'un raoddle statistique de 1' invention positicamg par le 
praticien dans une orientation approxiraative . La figure 3B reprg- 

25 sente le module 21' de la figure 3A 1' issue de l'6tape d' orien- 
tation sans modification de forme. La figure 3C repr^sente 
1' image 22 en trois dimensions resultant de la mise en oeuvre de 
1' invention, c'est-a-dire correspondant au modele 21', deforme 
pour que les rayons de rStrqprojection soient (idealement) tan- 

3 0 gents i tous ses contours. 

Par souci de simplification, un seul groupe de rayons r 
de retroprojection depuis une source 7 ont ete representes aux 
figures 3A a 3C et I'ecran de projection de 1' image correspon- 
dante n'a pas ete represent^. Comme il ressort des figures 3A a 

35 3C, un grand notribre de rayons r de la figure 3A ne sont pas tan- 
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gents aux ccxntours du moddle 21. Le nornbre de rayons de r6tropro- 
j action utilis6 depend de la precision souhaitee et des caracte- 
ristiques du modele statistique, essentiellement, du noctibre de 
points surfaciques choisis pour servir de reference lors de la 
5 definition de ce modele statistique. 

De la figure 3A a la figure 3B, la transformation est 
"rigide", c'est-a-dire que le modele n'est pas deforme. Lors de 
cette transformation, le modele 21 subit des translations et 
rotations afin d'obtenir des raesures d'ecarts minimum entre 

10 chaque rayon de retroprojection partant des points de contour de 
1' image et la surface de I'objet dans sa position courante. On 
utilise pour cela des methodes mathematiques de recherche du 
minimum de la sorame des carres de ces ecarts par rapport aux six 
parametres (trois degres de liberte en rotation et trois degr^s 

15 de liberty en translation) d^finissant la transformation rigide 
recherchie. On obtient vans configuration telle qu'illustr§e par 
la figure 3B dans laquelle le modele 21' a me position et 
orientation correcte mis una forme encore inparfaite. C'est 
pourquoi certain rayons r' ne sont pas tangents mais traversent 

20 le module. 

De pr6ference, la position resultant de la transforma- 
tion rigide (rotation, translation) est obtenue par une mgthode 
(algorithme) dite du point le plus proche (Iterative Closest 
Point, ICP) . Un exenple d'une telle methods est d^crit dans I'ar- 

25 tide "A Method for Registration of 3-D Shapes" de Paul J. Best 
et Neil D, McKay, paru dans IEEE TRANSACTIONS ON PATTERN ANALYSIS 
AND MACHINE INTELLIGENCE, vol. 14, NO. 2, en f^ier 1992, et 
dont le contenu est incoirpore ici par reference. 

De la figure 3B a la figure 3C, la transformation est 

30 "non- rigide" (elastique) et consiste a deformer le modele 21' 
sans changer son orientation jusqu'a obtenir 1' image 22 en trois 
dimensions pour laquelle tous les rayons de retroprojection r" 
sont a une distance rainimale des contours (idealement tangents) . 
Cette transformation non- rigide est effectuee, par exenple et 

35 selon un premier mode de realisation, en calculant les parametres 
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de d^fornation Bl partir d'vtn algorithms connu sous la denomina- 
tion angle- saxonne "Down Hill Sitrplex" . Get algorithme a ete 
decrit, par exarple, en 1965 par J. A. Nelder et R. Mead dans 
COMPUTER JOURNAL, vol. 7, pages 308-313, dont le contenu est 
5 incorpor§ ici par reference. 

Selon un mode de realisation prefere de 1' invention, 
1 ' algorithme utilise pour determiner la transfomation non-rigide 
du modele statistique est base sur la methode de Levenberg- 
Marquardt. Le principe de cette methode est de determiner le 

10 minimum d'une fonction non-lineaire multidimensionnelle en utili- 
sant les derivees partielles de la fonction par rapport aux para- 
metres de deformation du modele. Cette methode a ete decrite, par 
exemple, en 1963 par D.W. Marquardt dans JOURNAL OF THE SOCIETY 
FOR INDUSTRIAL AND APPLIED MATHEMATICS, vol. 11, pages 431-441, 

15 dont le contenu est incorpor^ ici par reference. La fonction dont 
on cherche le minimum est ici la somme des carr§s des distances 
entre un ensemble de rayons de rStrriprojection s'appuyant sur les 
points de contours et la surface du modele. Les pararaitres de la 
recherche du minimum sont les coefficients que I'on applique i 

2 0 chaque mode de deformation du module statistique, ainsi que cela 

est presents dans 1' article "Building a Cotrplete Surface Model 
from Sparse Data Using Statistical Shape Models: J^lication to 
Computer Assisted Knee Surgery" de Markus Fleute et Stgphane 
Lavallee, d^j^ cite. 
25 Cette methode est appliquee et repetee & chaque 

position approchee jusqu'a ce qu'il ne soit plus possible de 
minimiser les ecarts. On enchalne ainsi vn recalage rigide, un 
recalage elastique, un recalage rigide, un recalage elastique, 
etc. Une fois la forme en trois dimensions reconstituee, sa 

3 0 taille est bien entendu egalement connue dans la mesure ou toutes 

les vues en deux dimensions ont des tailles connues, de meme que 
la position, dans le referent iel de reference, du modele par rap- 
port aux sources. Le nombre de passes necessaire depend essen- 
tiellement du nombre de pcirametres cairacteristiques utilises pour 
35 modeliser la deformation. 
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On notera que I'ordre dans lequel sont effectu6es les 
deux §tapes ci-dessus est inportant. En effet, si on commence par 
effectuer vine deformation non-rigide, on va faire correspondre un 
models avec lane mauvaise orientation qu'il sera tres difficile a 
5 recuperer en raison de la deformation de forme subie. Toutefois, 
la minimisation des ecaits peut §tre effectuee globalement en 
prenant en cot^^te les parametres rigides et les parametres de 
deformations en meme tenps a chaque iteration du processus math6- 
raatique de minimisation. 

10 On notera que, i chaque etape iterative du precede dans 

laquelle on calcule les ecarts entre les contours du modele et 
les rayons de retroprojection, de nouveaux points entrant dans la 
mesure sont determines c6te modele et/ou cote rayons de projec- 
tion. Ceci const itue une distinction importante par rapport t des 

15 precedes connus de reconstitution d' images en trois dimensions 
dans lesquels les points de mesure sont les miSraes pour toutes les 
iterations . 

me autre caract^ristique de la pr^sente invention est 
que les Evaluations des hearts minimaux ne sont pas n^cessaire- 

20 ment faites de fagon exhaustive pour tous les elements de surface 
du module. En effet, cela conduirait parfois t des calculs trop 
dispendieux milme si le nottibre d' elements de surface peut etre 
minimise en utilisant des elements triangulaires . Selon 1' inven- 
tion, les evaluations d' ecarts sent faites par rapport t des 

25 lignes caracteristiques constituant des ggnSrateurs de contours 
du meddle. Cela minimise le nctribre de points de mesure. Par exem- 
ple, un OS de type femur, defini par environ 5000 bords de surfa- 
ces triangulaires, est defini par environ 300 generateurs de 
ccxitours . 

3 0 De preference, 1' analyse effectuee par le precede de 

1' invention ne touche que des generateurs de contour de 1' image 
en trois dimensions. Ainsi, la surface du modele statist ique est, 
de preference, formee d' elements triangulaires dont certaines 
aretes definissent des generateurs du contour. Un generateur de 

35 contour est defini par les aretes (bords) des triangles qui, pro- 
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jet^s sur vin plan, difinissent le contour (interne ou externa) . 
Le recours aux g6nirateurs de contours pentiet de reduire conside- 
rablement (par exettple, d'au moins un facteur 10) le ncrabre de 
points a chercher dans le inodele pour verifier la correspondance 
5 avec 1' image reconstituee. 

La figure 4 illustre la definition d'un generateur de 
contour d'une image (modele) en trois dimensions dent la surface 
est fomee d' elements de surface triaagulaires tels que decrits, 
par exettple, dans 1' article "Anatomy-based registration of ct- 

10 scan and intraoperative x-ray itnages for guiding a surgical 
rc±iot" de A. Gueziec, paru dans IEEE TRANSACTIONS ON MEDIO^ 
IMAGING, 17(5), pages 715-728, en octobre 1998, et dont le 
contenu est incorpore ici par reference. 

La figure 4 repr^sente, de fagon siitplifiee, deiix 

15 triangles 30 et 31 d^finissant une portion de surface d'un module 
en trois dimensions et dont I'argte cotnnune constitue une g€n§ra- 
trice 32 du contour. Le calcul t effectuer pour determiner si une 
argte constitue ou non une g^n^ratrice du contour consiste ^ cal- 
culer las angles a et P respectifs entre les normales aux surfa- 

20 ces des triangles 30 et 31 et les rayons 33 et 34 reliant le 
centre de ces surfaces au centre 35 de projection correspondant i 
la position de la source de rayons X. Si I'un des angles est 
inferieur t 90° tandis que 1' autre est supSrieur t 90°, leur 
arite cotnnune 32 est alors une g^n^ratrice du contour. 

25 La figure 5 illustre le type de tnesure effectu^e pour 

determiner I'^cart e entre un rayon r de retroproj action et un 
point d'un generateur de contour. Cette figure repr^sente, de 
fagon schittatique, une ittage I de I'objet dans un plan P et une 
forme 21 du modele statistique plac§e sur le trajet des rayons de 

3 0 rStrcprojection de 1' image I jusqu'i. la source (non representee) . 

Pour chaque rayon de retroproj action r, on cherche les 
points respectifs 41 et 40 du rayon de retroproj ect ion et d'un 
generateur de contour choisi parmi tous les gen^rateurs, pour 
lequel I'ecart e entre ces points est ttiiniitial. 
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Un avatitage de la prisente invention est qu'elle permet 
une reconstitution beaucoup plus rapide d'me image en trois 
dimensions par rapport aux techniques connues. 

Un autre avantage de la presente invention est qu'elle 
5 permet \m alignement correct du modele m§me dans des zones ou on 
ne dispose pas de donn^es de projection par les images en deux 
dimensions (par exeirple, certaines courbures internes) . 

On not era que le contour de surface pourra, le cas 
6cheant, etre affine au moyen d'un palpeur mecanique, optique ou 
10 magnetique foumissant des coordonnees spatiales dans le meme 
referentiel que les images en deux dimensions. On obtient alors 
des points supplementaires qui peuvent etre utilises dans la 
recherche de 1' image en trois dimensions. Le recours a un palpeur 
peut servir, par exeitple, a diminuer le nombre de vues en deux 
15 dimensions n^cessaires en donnant une information en trois dimen- 
sions . 

On notera ggaleraent que le precede de 1' invention peut 
s'appliquer ^ plusieurs contours surfaciques en trois dimensions 
constituant un ou plusieurs objets. Par exenople, on peut ainsi 

2 0 reconstituer siraultan§ment la surface exteme et la surface cor- 
ticale interne d'un os. On peut %aleraent reconstituer la forme 
de plusieurs os participant $l un articulation et visibles simul- 
taniment sur les images radiologiques (par exenple, le tibia et 
le f §mur si on s ' intSresse au genou) , ou encore la forme de plu- 

25 sieurs fragments ossevix d'lon mi§me os. Le module recherche doit, 
dans ce dernier cas, contenir les transformations rigides entre 
chaque contour surfacique en trois dimensions. Dans d' autre cas 
d' application i plusieurs c±)jets, il faut bien entendu connaitre 
les differents modeles statistiques et les Sventuelle transforraa- 

30 tions rigides ou elastiques entre eux. 

Bien entendu, la presente invention est susceptible de 
diver ses variantes et modifications qui apparaitront a I'homme de 
I'art. En particulier, le nombre de vues en deux dimensions 
devant etre utilise pour la reconstitution depend de la precision 

35 souhaitee et de la conplexite du modele. Dans certains cas, une 
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seule vue en deux ditnensicns pourra suffire. De plus, bien que 
1' invention ait ete dicrite en relation avec une source de rayons 
X mobile, celle-ci pourra etre rettplacee par plusieurs sources 
fixes, pourvu de respecter la contrainte d'obtenir des images en 
5 deux dimensions dans un metne rgfgrentiel. En outre, la mise en 
oeuvre de 1' invention fait bien entendu appal a des techniques de 
traitement numerique d' images dont la realisation pratique est a 
la portee de I'hctmie du metier a partir des indications fonction- 
nelles donnees ci-dessus et dans les pi±)lications cities en refe- 
1 0 rences . 
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REVENDICATIOMS 

1. Precede de reconstitution d'une image en trois 
dimensions representant les contours surfaciques d'au moins vn 
objet (1), a partir d'au moins une vue en deux dimensions de cet 
objet prise aux rayons X, caracterise en ce qu'il consiste §. : 

5 determiner la position de la source (7) de prise de 

vues dans vn r^f^rentiel de r6f6rence ; 

s^lectionner \jn module predefini constituant une forme 
moyenne (21) de 1' objet ; et 

de fagon it§rative jusqu'§. ce que les contours du 
10 module soient tels que les hearts, entre des rayons de r^tropro- 
jection des contours de 1' image en deux dimensions depuis la 
source et la surface du modele, soient minimaux : 

silectionner \ine orientation et xme position du modile 
dans le rifirentiel de r^f^rence, puis 
15 select ionner une deformation du modele pour modifier 

ses contours en trois dimensions. 

2. Proc4d§ selcn la revendication 1, caract4ris6 en ce 
que le module est dbtenu t partir d'une population d'cbjets pour 
laquelle on recherche la correspondance statistique commune ^ 

2 0 tous les objets pour determiner une forme moyenne et les deforma- 

tions principales par rapport k cette forme moyenne, de fagon ^ 
disposer d'au moins vm mcd§le statistique. 

3. Proc^dS selon la revendication 2, caract§rise en ce 
que les etapes de selection iteratives consistent a faire subir 

25 au modele statistique, successivement , une transformation rigide 
modifiant sa position et/ou son orientation et vine deformation 
non- rigide modifiant ses contours surfaciques. 

4. Procede selon I'lme quelconque des revendi cat ions 1 
a 3, caracterise en ce que les contours de 1' image en deux diinen- 

3 0 sions sont dbtenus automat iquement en projetant le module dans le 

plan de 1' image en deux dimensions, et en diformant les contours 
projetes de fagon i. les faire coincider avec les points de fort 
gradient en niveaux de gris de 1' image en deux dimensions. 
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5. Proc6d§ selon la revendication 4, caracterise en ce 
que la d^terraination autotnatique des contours de 1' image en deux 
dimensions est effectu^e de fagon iterative, chaque iteration 
etant intercalee entre devix iterations successives des etapes de 

5 selection. 

6. Precede selon I'une quelconque des revendi cat ions 1 
a 5, caracterise en ce qu'il consiste a determiner, dans le refe- 
rentiel de reference, des coordonnees en trois dimensions de 
points de I'objet, de fagon a disposer de points de reference 

10 supplementaires pour les etapes de selection iteratives de posi- 
tion, orientation, et deformation. 

7. Procede selon I'me quelconque des revendications 1 
Sl 6, caractSris^ en ce qu'il consiste a utiliser plusieurs images 
en deux dimensions pour lesquelles les positions respect ives de 

15 la source de prise de vues sent toutes determin^es dans le rgf§- 
rentiel de r§firence, et t effectuer les Stapes de s§lections 
itiratives en tenant compte des rayons de rStrcprojection des 
contours de toutes les images en deux dimensions. 

8. Proc^i selon la revendication 7, caract§ris€ en ce 
20 que le nottibre d' images utilise est fonction de la precision 

souhaitSe . 

9. Proc6d§ selon I'une quelconque des revendications 1 
k 8, caracterise en ce que la surface du modele est CQnstitu§e 
d' elements de triangle, lesdits ecarts etant mesures par rapport 

25 a des points de certaines ar6tes constituant des gen6rateurs du 
ccffitour en trois dimensions. 

10. Procede selon I'une quelconque des revendications 1 
S. 9, caracterise en ce qu'il est applique a la reconstitution des 
contours surfaciques de plusieurs objets lies entre eux par des 

30 relations de transformation rigide et/ou elastique. 

11. Procede selon I'une quelconque des revendications 1 
a 10, caracterise en ce qu'il est applique a la reconstitution 
d' images d'os. 



2798760 



21 

12. Systdme de traitement d' images, caract6rise en ce | 
qu'il coirporte des moyens pour la mise en oeuvre du precede selon 
I'Tone quelconque des revendications 1 11. | 
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